Insulator-Metal Transition in One Dimension 
Induced by Long-Range Electronic Interactions 
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The effects of a long range electronic potential on a one dimensional commensurate Charge Density 
Wave (CDW) state are investigated. Using numerical techniques it is shown that a transition to 
a metallic ground state is reached as the range of the electron-electron repulsion increases. In this 
metallic state, the optical conductivity exhibits a large Drude weight. Possible interpretations of 
our results are discussed. 

PACS numbers: 74.72.-h, 71.27.+a, 71.55.-i 



As the dimensionality of an electronic system de- 
creases, it is expected that charge screening would be- 
come less important in effectively reducing the range 
of the electron-electron interactions. Thus, in one di- 
mension (ID) the long-range (LR) Coulomb potential 
should play an important role in determining the physical 
properties of electronic models. 0] Experiments in GaAs 
quantum wires [|j and quasi-one-dimensional conductors 
H indeed highlight the importance of LR interactions 
between carriers. Recently, the role of a l/r Coulomb re- 
pulsion on the long distance properties of electrons con- 
fined to a chain has been theoretically investigated us- 
ing bosonisation techniques 0. The 4kp charge corre- 
lations decay very slowly with distance suggesting that 
the ground state (GS) is similar to a classical Wigncr 
crystal. However, this previous study was made in the 
continuum limit and the role of the lattice in this context 
is unclear. In addition, theoretical H and experimental 
H results suggest that Umklapp processes have a crucial 
importance in ID close to the Mott transition. 

In order to investigate the interplay between the 
Coulomb interaction and Umklapp scattering on the 
phase diagram and transport properties of ID systems, 
here static and dynamical observables are calculated on 
chains with a variety of electron-electron interactions us- 
ing numerical techniques. As the range of the electronic 
potential increases, an insulator to metal transition oc- 
curs. In particular, for unscreened l/r Coulomb inter- 
actions on a lattice our computational analysis suggests 
that the GS is metallic, an interesting result considering 
that for short-range interactions (SR) and the densities 
studied here the GS is a CDW insulator. 

The model analyzed here consists of a single chain with 
L-sites and an extended Hubbard-like interaction, 
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n)(iij — n), 
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where rij (n) is the local (average) electron density, and 
the rest of the notation is standard (t = 1 is the unit of 
energy). Vq — U/2 corresponds to the on-site interaction 
(U is the usual Hubbard local coupling) and will alway 
be included here. For \i — j\ > 1, three types of poten- 
tials will be mainly considered: (i) a plain LR Coulomb 
interaction V—j — Vc/\i — j\; (ii) a SR potential, typi- 
cally Vi-j = Vc/\i - j\ for \i - j\ < r max and Vi-j = 
otherwise, extending up to a distance r max (generally 
with r max = ro — 1 where ro is the average inter-particle 
distance, i.e. ro = 1/n in units of the lattice spacing); 
and (iii) a "mixed" potential with a Coulomb tail of ad- 
justable intensity, Vi-j = V c /\i - j\ for \i - j\ < r max 

■ max , < < 1. 
ID chains are 



and Vi-j = PV c /\i - j\ for \i - j\ > 



H = -*y)(ct a c i+ i jtT + h.c.) 



Using Exact Diagonalization techniques, 
here analyzed. A uniform ionic background has been in- 
cluded in Eq. (Q) such that the electrostatic energy per 
unit volume remains finite. 

Let us start with a purely ID commensurate CDW 
insulator (reference state) on which the LR Coulomb tail 
is switched on. For commensurate densities and with SR 
interactions, it is well-known that Umklapp processes can 
lead to an instability of the Luttinger Liquid GS towards 
an insulating state. Let us explore model Eq.(l) with a 
SR potential using r max = r$ — 1 = 2 at density n = 1/3 
varying the couplings U and Vc in order to determine a 
set of parameters for which the insulating nature of the 
GS can be unambiguously proven. Such insulating phase 
can be characterized by a vanishing charge stiffness: if the 
electrons experience a flux $ (in unit of the flux quantum) 
threaded through the ring, the GS energy Eq will become 
independent of $ once the chain length exceeds some 
characteristic localization length. In Fig. pl(a) the charge 
stiffness D — d 2 (E Q /L)/dK 2 (where k =^$) is shown 
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as a function of 1 /L 2 (a scalingbehavior also followed by 
Luttinger Liquid (LL) chains S). In the case of the SR 
potential the data strongly suggests an insulating state 
since D — ► when L — > oo for the couplings chosen here. 
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FIG. 1. Charge stiffness D vs 1/L 2 for n = 1/3. Dashed 
lines correspond to polynomial fits of the form a + b-^ + Cj^ . 
(a) Comparison between SR and LR potentials for open shell 
configurations. Lhe SR potential extends up to r max — 2; (b) 
Comparison between Vc m (r) (filled symbols) and 1/r (open 
symbols) potentials for boundary conditions giving rise to 
open or closed shell configurations. 

The insulating character of this state can be fur- 
ther established from the finite size scaling behavior of 
the single particle excitation energy (Fig.g), defined as 
A c ,i( L ) = E o(nL + 1, L) + E (nL -l,L)- 2E Q {nL, L), 
where Eq(N £ , L) is the GS energy for N e electrons and L 
sites H . For the SR potential and the same parameters as 
previously, Fig. ||(a) strongly suggests that Ac,i extrap- 
olates to a finite and sizable gap. In addition, explicitly 
evaluating real space charge correlations in this state a 
well-pronounced 4fci? = 2nh pattern was observed. The 
physical nature of this state is then clear: it forms a com- 
mensurate CDW state in which the particles are localized 
on every ro sites. In the thermodynamic limit, the GS is 
r -fold degenerate with gapped charge excitations pL| ]. 

Let us now consider a more extended interaction. The 
U = oo phase diagram of the repulsive V1—V2 model at 
quarter filling (n — 1/2) obtained from the investigation 
of the charge gap Ac t i is shown in Fig. || (V\ and V2 are 
the electronic density-density repulsion at distance 1 and 
2, respectively, while at larger distances the interaction 
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FIG. 2. Finite size scaling of the single particle gap Ac,i 
using U = 16 and Vc = 10. Dashed lines correspond to 
polynomial fits of the form a + bj- + Cj^. (a) Comparison 
between SR and Coulomb (LR) interactions, (b) 1/r Coulomb 
potentials. The positions of the absorption peaks labeled by 
I and II in Fig. 4(b) are also shown (open symbols). 



is zero). For such a density (ro = 2), the line V% = 
corresponds to an insulating CDW state above a criti- 
cal value of V\ in agreement with our previous results 
for n — 1/3, and with earlier calculations using spinlcss 
fermions by Emery and Noguera jllj . When a finite V2 is 
slowly switched on (now r max — ro) the insulating phase 
is destroyed above a critical line. Thus, a repulsion be- 
yond distance ro — 1 drives the system metallic starting 
from a CDW insulator. Of special importance is the case 
V2 = Vi/2, which fulfills the relation between the first 
terms of the Coulomb interaction, since it corresponds to 
a particularly stable metallic GS. It is also important to 
stress that the metallic state is here of the Luttinger Liq- 
uid (LL) type with uj p (q) ~ u p q charge excitations. The 
corresponding values of the LL exponent a characterizing 
the density of states || (N(uj) ~ u> a ) are shown on Fig. 3. 
This exponent was obtained using a = (K p + K~ 1 — 2) /4, 
and K p was calculated from the compressibility and the 
charge velocity following a standard procedure. Q 

Let us now turn to the case of a 1/r Coulomb potential. 
To facilitate the comparison with the previous results 
for the CDW insulating regime, the density n = 1/3 at 
U = 16 and Vc = 10 will be considered. When the 
LR tail is added, in contrast to the SR potential, the 
scaling of the charge stiffness in Fig. [l](a) suggests that it 
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FIG. 3. Phase diagram of the V1-V2 model (U = 00) at 
a density of n = 1/2. See also Ref. jllj]. The extrapolated 
values of the exponent a are indicated on the plot. 



reaches a finite value as L — ► 00. To check the accuracy 
of the finite size extrapolation the sine approximation to 
the Coulomb potential V^ m (r) = ^t/ sin (r/nL) was also 
considered. A comparison between Vq" and the plain 
1/r interaction is shown in Fig. |l|(b). All data sets are 
consistent with the same bulk extrapolation. Then, the 
present results show that the commensurate CDW state 
undergoes a transition to a metallic state as the range of 
the electronic interaction grows. 

The metallic character of the 1 Jr Coulomb model can 
also be further established from the finite size scaling of 
the single particle gap Ac,i as shown in Figs. ||(a,b). For 
the sets of parameters considered here, a robust 1/L scal- 
ing behavior is found which suggests that the charge gap 
vanishes in the bulk, consistent with the finite value of 
the charge stiffness found before. It is interesting to ob- 
serve that using a weaker Coulomb tail of amplitude (3V C 
(/3 < 1) also drives the system into a metallic state as 
shown for (3 = 0.5 in Figs. |](a) and ||(a). This suggests 
that the range of the interaction, rather than its inten- 
sity, is the key element in inducing the insulator-metal 
transition. Also note that the model with LR interac- 
tions is metallic apparently at any value of Vc/t, i.e. an 
infinitesimal amount of kinetic energy is enough to de- 
stroy the insulating state. This is to be contrasted with 
SR models, such as the t — V\ model, where an insulator- 
metal transition occurs only for a large enough hopping 
amplitude |pj| . In addition, the metallic state exhibits 
collective spin modes similar to those of a LL, i.e. with 
a dispersion <jJ a {q) ~ u a q. Using finite size chains u a 
can be easily calculated on closed shell configurations for 
L = 12, 18 and 24 sites, by means of the energy differ- 
ence between the first triplet state with momentum 2tt/L 
and the GS. We have found that u a follows very closely a 



1/L 2 scaling behavior, which enables an accurate deter- 
mination of its bulk value (table I). The small values of 
u„ indicates that the effective exchange coupling between 
neighboring particles at a distance rn is small. 

To gain insights on the transport properties of the 
metallic state observed here, the zero-temperature op- 
tical conductivity 



<j(u>, T = 0) = 7rD 6(ui) + a re g(oj), 



(2) 



was numerically investigated. Results for the LR poten- 
tial are shown in Fig. [|(a) for L = 18 and L = 24 rings. 
The shape, as well as the magnitude of the absorption 
curve, depends weakly on the system size. By calculat- 
ing separately a reg (ui), and using the previous results for 
the charge stiffness D, it can be explicitly checked the 
validity of the optical sum rule nD/2 + I = ^E^/L, 
where E^/L is the kinetic energy per site and Iq is the 
finite frequency integrated absorption. Extrapolations 
(L — > 00) of these quantities are provided in Table I. 
It is clear that the metallic and insulating regimes have 
very different transport properties. Specially, the weight 
of the absorption band of the metallic state (given by 
<T reg (u>)) corresponds only to 1 to 2% of the total weight, 
while it exhausts the total sum rule in the case of the in- 
sulator |D| . On the other hand, the energy scales of the 
absorption bands are comparable for the insulator and 
the metal (Fig. |(b)). With increasing system size the 
largest low energy peak in a reg (uj) shifts to lower energy 
(see a comparison with Ac,i in Fig. ||) as its weight de- 
creases. These excitations might correspond to excitons 
which have lower energies than the CDW gap when finite 
range interactions are considered. 

A simple interpretation of the insulator-metal transi- 
tion observed here is as follows: consider the potential 
experienced by a mobile electron on a chain assuming the 
rest of the electrons are fixed into a CDW configuration, 
with particles equally spaced at distances rn = 1/n. For 
a strong SR interaction the mobile electron is trapped in 
a deep square-well potential, compatible with the insu- 
lating properties observed numerically However, as the 
potential range increases the potential experienced by the 
mobile electron diminishes. Actually, in the (unphysical) 
limit where the repulsive potential between electrons is 
made distance-independent such potential becomes irrel- 
evant. This "one-electron" interpretation is favored by 
the fact that the insulator-metal transition is not associ- 
ated only with the (usually subtle) 1/r interactions but 
it appears also for the SR V\ — V2 model as well (Fig. ||). 

An alternative interpretation of our results involves a 
transition from a pinned CDW to a Wigner crystal driven 
by the LR part of the potential juj . The bosonisation 
calculations suggest that the low frequency, long wave- 
length conductivity of the Wigner crystal should behave 



<j(u),q) oc 



i(u> + ie) 



(u + ie) 2 - {u p {q)f 



(3) 
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excitations in this scenario. 

Quasi-ID organic conductors @ like (TMTSF) 2 PF 6 
are compounds where 1 jr Coulomb repulsion mi ght play 
an important role. This is suggested by ab-initio [|l7| and 
Hartree Fock- Valence bond |L8| quantum chemical cal- 
culations. Note also that despite its many successes, the 
description of quasi-one-dimensional conductors based on 
the LL approach still faces various problems to provide a 
consistent explanation of all the experimental data such 
as NMR relaxation rates JhJ, photoemission || and IR 
absorption |6j. As argued in this paper, a long-distance 
1/r Coulomb repulsion within the chains is expected to 
have drastic effects on the transport properties of such 
ID systems, and this phenomenon may contribute to a 
more accurate theoretical description of ID conductors. 
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FIG. 4. Optical absorption of a 1/r-Hubbard chain, (a) 
Data for 18- and 24-sites chains, (b) Comparison with the 
insulating CDW state (/3 = 0). The small absorption in the 
case of the LR Coulomb potential (/3 = 1) has been multiplied 
by a factor 50 and the comparatively (very) large Drude peak 
at uj = has been omitted for clarity. 

In general, the q — > and lu — > limits do not com- 
mute. The relations cr(uj,q = 0) = 2Z?Drudc<5(w) and 
<j(0,q) = G5(q) define the Drude weight DDrude and the 
conductance G, respectively. Using uj p {q) oc g| In g| 1,/2 
(plasmons), Q Eq. (|J) leads to a finite Drude weight and 
a vanishing conductance. The Drude weight is directly 
related to the charge stiffness calculated previously by 
-DDrude = 7T-D/2. Therefore, our numerical results for the 
LR potential on the lattice are also compatible with the 
bosonisation approach in the continuum limit, suggest- 
ing that a metallic Wigner crystal GS can be realized 
on a lattice. |l5| The metallicity is caused by collective 
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u = 


10, v c = 


5/3 = 1 


0.13(5) 


0.836 


0.010(3) 


0.846 


u = 


16, V c = 


10 p = 1 


0.04(1) 


0.752 


0.017(1) 


0.769 


u = 


16, V c = 


10 p = 


0.00(0) 


0.0 


0.43 


0.43 



TABLE I. Extrapolated values of the spin velocity u a , the 
Drude weight, the integrated optical absorption and the total 
sum rule for n = 1/3. 
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